Quantum master equation for a system influencing its environment 
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We derive a new perturbative quantum master equation for the reduced density matrix of a system 
interacting with an environment (with a dense spectrum of energy levels). The total system energy 
(system plus environment) is constant and finite. This equation takes into account the finite energy 
effects of the environment due to the total energy conservation. This equation is more general than 
the common perturbative equations used for describing a system in interaction with an environment 
(like the Redfield equation 0, I3 or the Cohen- Tannoudji one Q) because these last equations 
can be deduced from it in the limit of an infinitely large environment. We apply numerically this 
equation to the spin-GORM model. This model represents the interaction of a two-level system 
with an environment described by random matrices. We compare our equation with the exact von 
Neumann equation of the total system and show its superiority compared to the Redfield equation 
(in the Markovian and non-Markovian cases). 
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The study of the quantum dynamics of a system inter- 
acting with its environment is an old but still very impor- 
tant problem in quantum mechanics and more specifically 
in non-equilibrium quantum statistical mechanics. The 
understanding of such problems is central for the study 
of very fundamental physical processes such as the relax- 
ation to the thermodynamic equilibrium, the decoherence 
or the thermalization at the environment temperature. 
Actually, one of the more general and of the most used 
equations to describe the dynamics of a system interact- 
ing with its environment is the non-Markovian Redfield 
equation. The other well-known equations used in the 
past for the study of such systems (the Markovian Red- 
field equation 1], the Cohen- Tannoudji equation or 
the Lindblad equation A\) can be derived from the non- 
Markovian Redfield equation 0. The Markovian form 
of this equation was first derived empirically for under- 
standing NMR experiments but later generalized to 
the non-Markovian case by deducing it formally by per- 
turbation theory from the von Neumann equation of the 
total system (system plus environment) The funda- 
mental assumption in the Redfield theory is that the envi- 
ronment is infinite and therefore not affected by the sys- 
tem. As a consequence of the recent development of nan- 
otechnology, there is an increasing interest for nanometric 
systems. In nanosystems, one can encounter situations 
where purely quantum degrees of freedom like the spin 
are interacting with degrees of freedom that are dense 
in energy but with an energy distribution that varies on 
small energy scales sometime of the order of the system 
energy. In such situations, the Redfield equation is not 
valid anymore, the infinite environment hypothesis fails 
and the finite energy effects of the total system have to be 
taken into account. In this paper the plan is as follows: 
we derive our new equation, that takes into account the 
finite energy effects of the total system, we then apply it 
to the population dynamics of the spin-GORM model (a 
two-level system interacting in a non-diagonal way with 



an environment described by Gaussian orthogonal ran- 
dom matrices) and we finally compare it to the exact 
von Neumann equation for the total system and to the 
Redfield equation. We will finally conclude presenting 
the perspectives created by this work. The quantum sys- 
tem that we consider has a simple spectrum (integrable 
system) given by his Hamiltonian Hs- Its eigenvalues 
(respectively its eigenvectors) are given by {Eg} (respec- 
tively by {|s)}). The environment is also a quantum sys- 
tem but with a dense spectrum (containing many levels 
forming a quasi-continuum) given by a Hamiltonian Hb ■ 
Its eigenvalues (respectively its eigenvectors) are given 
by {e} (respectively by {|e)}). The interaction between 
the system and the environment is generally taken to 
be the product of a system operator S and an environ- 
ment operator B. The coupling parameter A measures 
the intensity of the interaction between the system and 
the environment. Therefore, the Hamiltonian of the to- 
tal system is given by i?tot = Hs + Hb + XSB. The 
exact time evolution of the total system is described by 
the von Neumann equation p{t) — —i[Htot, pit)], where 
p{t) is the density matrix of the total system. The sys- 
tem dynamics is described by the reduced density matrix 
of the system ps{t) = Ti Bp{t)- The total system has a 
finite constant energy. At initial time, the environment 
is in a microcanonical state at energy e. 

The ansatz in the Redfield derivation is to suppose that 
the total density matrix evolves keeping the following 
form: 



p{t)^ ps{t)^6{e-HB). 



(1) 



We see that the environment part of the density matrix 
does not evolves, supposing that the environment is not 
affected by the dynamics. But as we announced it, we 
want to include the finite total energy effects in the dy- 
namics. The main idea (the new ansatz) is to suppose 
that the total density matrix can be described at all times 
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by a density matrix of the following form: 
n{HB) 

where n{e) = TrB(5(e — Hb) is the smoothed density of 
states of the environment. The reduced density matrix 
of the system becomes ps{t) = J deTis^i^ — HB)p{t) — 
Tlis s' / dePss'i^', t). The only approximation made 

by the ansatz is that we neglect the contribution to the 
dynamics coming from the environment coherences. We 
see that the environment energy can now depend on the 
system state. Inserting ^ in the von Neumann equa- 
tion, taking the trace over the environment degrees of 
freedom and performing a perturbative expansion up to 
the second order in A, one gets our new equation. For 
the population dynamics, this equation takes the follow- 
ing form: 



s,s' 

+{s\S\s'){s'\S\-s) 
Pss{£;t) / de F{e,€ )n[e )- 



(3) 



e - e' 



-{s\S\s){s'\S\s) 
n(e) / de F{e,e)Pss'[e ;t) 



Es - E-s, + e - e' 

where F{e,e') = "|(e|S|e')p' where the quotes denote 
a smoothening over the dense sppectrum of eigenvalues 
around e and e' . Performing the same procedure using 
instead of ^ gives the non-Markovian Redfield equa- 
tion ^ which can be seen as the particular case of our 
equation when the environment density of states varies 
on a large energy scale compared to the typical energies 
of the system. The Markovian approximation consists in 
taking the infinite time limit of the time-dependent coef- 
ficients of lj2Jl using the property limi-^oo = ttS{(_). 
Performing this approximation on our new equation (re- 
spectively on the non-Markovian Redfield equation) gives 
the Markovian version of our new equation (respectively 
the Redfield equation ^^l). If one further neglects the con- 
tributions of the coherences to the populations evolution 
and the contributions of the populations to the coher- 
ences evolution, one gets a simplified Markovian version 
of our new equation. This equation takes the following 
form for the populations dynamics: 

Pss{e;t) = -27rA2^ |(s|^|s')pF(e,S,-i;,- -l-e) 

n{Es- Es' +e)Ps,s{e;t) 

n{e)Ps,,s'{Es- Es, +e;t){i) 

Doing the same for the Redfield equation one get the 
well known Cohen- Tannoudji equation Our simplified 



Markovian equation describes the evolution of the total 
system as a random walk between the states belonging to 
the same energy shell with transition probabilities given 
by the Fermi golden rule. This equation may look like 
the Pauli equation 6] but it describes the time evolution 
of the distributions of populations over the energy of the 
environment which is the new feature of our equation. 
On the other hand, the difference between our simpli- 
fied Markovian equation (respectively the non-Markovian 
version of our equation) and the Cohen- Tannoudji equa- 
tion (respectively the non-Markovian Redfield equa- 
tion) is the fact that this last equation considers that 
the environment density of states is not affected by the 
system energy. This is represented in Fig. ^ We will 
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FIG. 1: Representation of the energy exchanges described 
respectively by the Cohen- Tannoudji and by the simplified 
Markovian version of our new equation for a four level sys- 
tem. The system energy is represented on the abscissa. The 
environment energy (continuum) is represented on the ordi- 
nate. The width in energy of the environment spectrum is 
fixed to 1 and its density of states is supposed semicircular. 
The total energy of the system is given by E. The initial con- 
dition is denoted by two empty superposed circles. We see 
that transitions that preserve the energy of the total system 
have to occur along the total energy line crossing the plane. 
Doing this, they satisfy the Fermi golden rule for the total 
system. One can see that only our equation satisfy this con- 
dition. The Cohen- Tannoudji equation describes transitions 
that occur along a vertical line at constant environment en- 
ergy and is therefore wrong when the system energies are of 
the order or larger than the typical energy scale of variation 
of the environment density of states. 

apply now our equation to the spin-GORM model. This 
model describes the evolution of a two-level system that 
interacts in a non-diagonal way with a complex environ- 
ment. Here, the complexity is supposed to come from 
many-body interactions like in heavy nuclei or from a 
classically chaotic dynamics. Therefore, we represent 
all the environment operators by Gaussian orthogonal 
random matrices (GORM). The Hamiltonian of the to- 



tal system is then Hu 
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where X and X' are Gaussian orthogonal random ma- 
trices of size N and probabiUty density proportional to 
exp{—jTrX^ . The three fundamental parameters deter- 
mining the model are A, A, and N. The smoothed den- 
sity of states of the environment is given by the Wigner 

semicircular n(e) = ^ — with the convention that 

y/x = for a; < 0. The smoothed density of states of the 
nonperturbed total system (A = 0) is the sum of two 



and 



A 



semicircular centered at the two energies - 
depicted in Fig. |21 
In this paper, we restrict ourselves to the study of the 
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FIG. 2: Schematic representation of the smoothed density of 
states of the nonperturbed total system (A = 0) for different 
values of A. 



the lower one), one gets Q 

/&'[P++(e';t)-P__(6' + A;0] (5) 



where 



(6) 



--P±±{e;t) 



A^ 

TT 



de'P^^{e'-t) 



(±A + e - e') 
sin(±A + e - e')t 



(±A- 



An important remark is that the Markovian and the sim- 
plified Markovian version of our equation are the same for 
the spin-GORM model. Neglecting the coherences con- 
tributions to the populations dynamics is not necessary 
here. This has also the consequence that the Marko- 
vian Redfield equation and the Cohen- Tannoudji arc the 
same. 

The Markovian version of our equation 10) becomes 
the following equation for the spin-GORM model: 
{a^)^{t) = P++{e;t)-P^^{€ + A]t). In this case, we get 
an analytical solution describing an exponential decay to 
equilibrium: 



where the relaxation rate is given by 



"ft , 



(7) 



small coupling regimes A <C 1 because all the equations 
discussed in this paper are obtained perturbatively. An- 
other restriction is related to the mean level spacing of 
the environment. The coupling between the levels (that 
is of order A^ because the first order in perturbation the- 
ory is zero due to the non-diagonal nature of the cou- 
pling) has to be of the order or larger than the mean 
level spacing of the environment to induce a sufficient 
interaction between the levels belonging to the micro- 
canonical energy shell in order to reach a microcanonical 
distribution inside these shells. The criterion is therefore 
> Of course, the lower bound disappears in the 
continuum limit N oo. The validity domain is shown 
in Fig. 13 The detailed study of the lower bound of the 
coupling parameter has been done in . 
Applying our equation Q to the spin-GORM model in 
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FIG. 3: Schematic representation of the validity domain of 
the kinetic equation in the coupling parameter A. 

order to study the system population evolution through 
tJj, (difference between the probability of beiirg in the up- 
per state of the system minus the probability of being in 



(e + A) 



and the equilibrium population by 



7 )^ 



{e + Ay 



Notice that the relaxation rate and the equilibrium pop- 
ulation are independent of the initial condition. 
We now start the discussion of the comparison between 
the different equations based on the numerical simula- 
tions of Fig. 21 We always take ((Tz)^(O) = 1 in our 
numerical simulations. The initial state of the environ- 
ment is characterized by the energy of the microcanonical 
distribution e. The width of the energy shell is always 
Se — 0.05. The curves are averages over % = 10 realiza- 
tions of GORM of size N = 2000 of the Hamiltoniair, and 
over the different eigenstates of the environment energy 
shell. 

As we announced it, our equation is necessary when 
the system energy A is of the order of magnitude of the 
energy scale of variation of typical environment smooth 
density of states: n{e + A) ^ '^(c)- But when A is small 
enough and therefore n(e + A) w '^(c)? our equation re- 
duces to the non-Markovian Redfield equation. This can 
be well seen comparing Figs. and B. In Fig. 0J\, our 
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FIG. 4: Two-level system dynamics for the spin-GORM 
model. 



Markovian (M) and non-Markovian (NM) equations are 
very close to the Redfield M and NM ones, but this is 
not the case anymore in Fig. 0j3. Another important 
point is the validity of the Markovian approximation. 
This approximation has the effect of constraining the dy- 
namics inside the microcanonical energy shell of the non- 
perturbed spectrum according to the Fermi golden rule. 

At short time, because ^'"^A+e-l'V* ^ delta 

function, the non-Markovian equation has the possibility 
of describing a spread of the probability in energy around 
the microcanonical energy shell. When the delta contri- 
bution to the dynamics is important, the non-Markovian 
effect are very small, as in Figs. and B, and can 
only be seen on very short time scales. But when this 
contribution is small or zero, like in Fig. 2P (because 
the microcanonical energy shell is not inside the overlap- 
ping zone as seen in Fig. [2J), then the non-Markovian ef- 
fects become very important. In Fig. ^P, the Markovian 
curves are not represented because they completely miss 
the dynamics (they predict no evolution (a^) — 1). Fig. 
0P represents therefore a pure non-Markovian dynamics 
to which only the non-central part of '^'"jlA+e-'e'V* 
tributes. 

We can conclude saying that our new equation is always 
valid in the small coupling limit independently of the en- 
ergy ratio between the system and the environment. It is 
therefore an important equation for the study of nanosys- 
tems. This equation reduces to the Redfield equation 
for very small system energies compared to the energy 
variations of the environment density of states. It can 
be shown that it is also in this limit that the equilib- 
rium values of the system populations thermalize to a 
canonical distribution corresponding to the microcanon- 
ical temperature of the environment 0. The coherences 
dynamics will be investigated in future work. 
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